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Abstract: DNA sequences from the mitochondrial gene cytochrome oxidase subunit I (mtDNA CO 7) were used to 
estimate the genetic variability in two wild populations and two cultured populations of the hard shelled mussel, Mytilus 


coruscus. Thirty haplotypes were identified in the four populations. The cultured populations exhibited a lower number of 


haplotypes and genetic diversity than those of the wild populations, suggesting that a small number of effective founding 


breeders contributed to the genetic variation of the cultured populations. No significant differentiation was observed 


between the cultured population and local wild population, implying that persistent gene flow occurred in these 


populations. This genetic survey is intended as a baseline for future genetic monitoring of M. coruscus aquaculture stocks. 
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厚 壳 贻 贝 养殖 群体 与 野生 群体 线粒体 DNA 的 遗传 分 析 
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(上 海 海洋 大 学 水 产 与 生命 学 院 ， 上 海 201306) 
摘要 : 采用 线粒体 DNA COI 基因 序列 对 厚 壳 贻 贝 2 个 养殖 群体 与 2 个 野生 群体 的 遗传 多 样 性 进行 了 研究 。 
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个 群体 共 获 得 30 个 单 倍 型 。 结 果 显 本 
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: 在 养殖 群体 中 单 倍 型 的 数量 和 遗传 多 样 性 要 比 野 生 群 体 的 低 ， 可 能 是 
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只 有 少量 的 有 效 父母 本 对 养殖 群体 的 遗传 变异 有 贡献 所 致 。 养 殖 群体 与 当地 野生 群体 之 间 也 未 发 生 显著 的 遗传 
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化 ， 可 能 是 因为 它们 之 间 存 在 基 
从 而 保证 养殖 群体 的 遗传 多 样 性 水 平 。 
关键 词 : 厚 壳 贻 贝 ; 
中 图 分 类 号 : Q959.215; Q349 


> 






































文献 标识 码 : A 


The hard shelled mussel, Mytilus coruscus, is an 


economically important mussel which is widely 
distributed from coast of China, Japan, and Korea (Wang, 
1997). It is documented that this species occurred from 
Dalian to Xiamen in China. The culture of M.coruscus 
has been carried out in only a few regions of China in the 
past decades. The cultured juveniles mainly originated 
from the collection of natural populations. In recent years, 
natural juveniles have decreased while mussel farmers 
have increased. Due to overexploitation, most mussel 
stocks have dramatically declined. Therefore, it is very 
important to obtain juveniles that are produced by 


artificial breeding to increase the juvenile population. 
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流 。 在 今后 厚 壳 贻 贝 养殖 过 程 中 ， 本 研究 可 以 用 在 对 养殖 群体 进行 遗传 监测 


群体 分 化 : 遗传 多 样 性 ， 线 粒 体 DNA COI 基因 






























































文章 编号 : 0254-5853(2009)03-0240-07 


The first hatchery stock is being developed in Zhejiang 
China. 
frequencies between farmed strains and the wild source 


province, However, differences in allele 
populations have been shown to be a result of breeding 
related individuals or the use of small numbers of 
individuals as brood-stock, leading to lower genetic 
variability of the farmed strains (McGinnity et al, 2003). 
Additionally, farmed strains of small population size are 
more sensitive to genetic drift and consequently have 
lower within population genetic diversity than wild 
populations (Allendorf, 1986). Genetic variability is the 
primary resource in the successful artificial propagation 


of any species. Proper management and breeding 
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programs must be implemented to preserve genetic 
variability, and prevent inbreeding depression. However, 
for such programs to be successful, information on the 
genetic relationships among cultured and wild 
populations is essential. 

Mitochondrial DNA has been used extensively for 
2000). In 


comparison to allozyme or nuclear DNA, the higher 


studies in molecular ecology (Avise, 
mutation rate, smaller effective population size is 
expected to provide greater power to identify population 
structure. Additionally, mtDNA gives a better estimate of 
genetic differentiation than nuclear markers since it is 
approximately fourfold more sensitive to genetic drift 
and founder effects (Birky et al, 1983). Mitochondrial 
Cytochrome Oxidase I gene (mtCOI) has been used 
widely in marine species, such as Tegillarca granosa 
(Zheng et al, 2009), due to adequate levels of variability 
and easy amplification via universal primers (Folmer et 
al, 1994). 

In this study, we used the mitochondrial 
Cytochrome Oxidase I gene (mtCOI) to describe the 
genetic variability of two cultured, and two wild 


populations and also to quantify the genetic differentiat- 
ion between them. 


1 Materials and Methods 


1.1 Samples 

Specimens of M.coruscus obtained from both 
cultured and natural populations are listed in Tab. 1 and 
Fig. 1. The two wild populations were randomly sampled 
by divers from the subtidal zone of Shengsi and 
Lianjiang in 2007, respectively (coded WSS and WLJ). 
The two cultured populations were of Shengshan (coded 
CSS) and Huaniao (coded CHN), and also collected in 
2007. The CSS population was collected from the 
cultured area of M.coruscus. In this area, the aquaculture 
practice had been carried out for three decades. A large 
number of the natural spat collected were commonly 
reared. The CHN population hatchery stock was the 
first-generation of offspring produced in spring 2005 
using hundreds of wild caught M.coruscus in Shengsi 
(exact numbers not documented). Gill tissues were 
obtained and stored in 100% ethanol at room temperature 
until DNA extraction. 


Tab.1 Wild and cultivated mussels (Mytilus coruscus) obtained at different sites 


used in this study 





Origin Site Abbreviation Longitude Latitude Year 
Wild Shengsi WSS 122? 41" 30* 51" 2007-11 
Wild Lianjiang WLJ 119° 53’ 26°02’ 2007-05 
Cultured Huaniao CHN 122? 40' 30° 52' 2007-05 
Cultured Shengshan CSS 122? 48' 30° 43' 2007-05 


12 Extraction of DNA 

Genomic DNA was isolated from gill using the 
standard proteinase K digestion and phenol/chloroform 
extraction procedures described by Shen et al (2006). 
DNA quality was assessed by running samples on 1% 
agarose gels, and DNA concentration was measured with 
an UV/visible spectrophotometer (Eppendorf AG 22331 
Hamburg) for absorption at 260nm. DNA was diluted to 
50 ng/uL for reaction (PCR) 
amplifications. 


polymerase chain 
13 Amplification of DNA and sequencing 

PCR was carried out in 50uL reactions containing 
100ng DNA sample, 5uL 10xPCR buffer, 2.0 jumol/L 
MgCb, 200 umol/L dNTPs, 0.2 umol/L primer, 1U Taq 
polymerase. A pair of universal primers was used in this 
study (forward primer LCO1490: 5'-GGTCAACAA- 
ATCATAAAGATTGG-3'; reversal primer HCO2198: 
5'-TAAACTTCAGGGTGACCAAAAAATCA-3'(Folmer 
et al, 1994). PCR was performed using the Eppendorf 


PCR system programmed for an initial denaturation step 
of 3min at 94°C followed by 35 cycles, each consisting 
of 94°C denaturation for 50 second, 48°C annealing for 
50 second and 72°C extension for 1 min. A final 
extension of 10min was performed at 72°C and the PCR 
products were then held indefinitely at 4°C. A negative 
control, consisting of all the reaction components except 
template DNA, was also included for each of 
amplification. 

PCR products were visualized using 1.596 agarose 
gel stained with ethidium bromide. All amplified 
products were purified using TIANquick Midi 
Purification Kit (Tiangen, China). Purified PCR products 
were directly sequenced in both directions using the PCR 
primers on an Applied Biosystems ABI 3730 DNA 
sequencer. 

1.4 Data analysis 

For all sequence analyses, sequences were aligned 

with BioEdit Sequence Alignment Editor version 7.0.9 
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Fig. 1 Sample locations of two wild and two cultured Mytilus corcuscus 


(Hall, 1999) and saved in the Fasta format. The identical 
haplotypes in the aligned matrix were identified and 
collapsed using COLLAPSE version 1.2 (Posada, 2006). 
Nucleotide composition was computed in MEGA version 
4.0 (Tamura et al, 2007). Genetic diversity among wild 
and cultivated populations was estimated using 
haplotype number (N), gene diversity (A), and nucleotide 
diversity (x) as implemented in DnaSP version 4.0 
(Rozas et al, 2003). Gene diversity was calculated 
according to Nei (1987), using the probability that two 
randomly chosen haplotypes of the sample are different. 
Nucleotide diversity was also computed according to Nei 


(1987), as the probability that two randomly chosen 
homologous sites are different. Genetic diversity was 
compared between farmed and wild populations of 
mussel using FSTAT (Goudet, 1995). 

Pairwise estimates of Fs and corrected average 
pairwise population distances among all populations 
were obtained using the program Arlequin3.1 (Excoffier 
et al, 2005). The significance of these estimates was 
tested by comparing observed Fx with a null distribution 
obtained by 10000 random permutations of the data set 
(Excoffier et al, 1992). 


2 Results 


No. 3 


2.4 Characteristics of mtDNA COI of M.coruscus 

Once all sequences were aligned, a total segment 
length of 573bp of cytochrome oxidase subunit I (COT) 
was obtained. The average nucleotide composition in our 
samples is A=33.6%, C=23.9%, G=15.8%, T=26.6%, 
with an average G+C content of 39.7%. These sequences 
include 59 variable sites, of which 41 are singleton 
variable sites and 18 are parsimony informative sites. 
Sequences have been deposited in GenBank under the 
following accession number range: FJ495257-FJ495286. 
2.2 Genetic variability and population structure 

within populations 

The samples examined show a wide range of values 
for the number of haplotypes, gene diversities and 
nucleotide diversities (Tab. 2). A total of 30 haplotypes 
were obtained for wild and cultured populations of M. 
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coruscus (Tab. 3). The number of individuals of different 
haplotypes ranged from 1 to 18. The haplotype Hap2 had 
the greatest number of individuals across all samples. 
The wild WSS population had the greatest number of 
haplotypes, whereas the cultured CHN population 
originating from the Shengsi area had the fewest number 
of haplotypes. Every population had population-specific 
haplotypes. In the cultured WLJ population, only 3 
specific haplotypes were found, exhibiting a decrease of 
specific haplotypes. In contrast, the samples from the 
wild WSS population exhibited the greatest number of 
population-specific haplotypes. The samples from natural 
populations exhibited the greatest number of haplotypes 
and population-specific haplotypes compared to cultured 
samples. 


Tab. 2 Genetic diversity parameters of Mytilus coruscus populations examined 





Population Number of individuals(n) Number of haplotypes(N) Gene diversity(/) Nucleotide diversity(z) 
WSS 20 13 0.911 0.00898 
CHN 13 6 0.718 0.00293 
CSS 20 12 0.911 0.00834 
WLJ 13 10 0.923 0.00948 
Overall 66 30 0.895 0.00832 


The samples from wild WSS population had an H 
value of 0.911 and z of 0.00898, compared to the 
aquaculture samples CHN and CSS with H values of 
0.718, 0.911 and z values 0.00293, 0.00834, respectively. 
The wild WSS population had 3-fold higher nucleotide 
diversity than the cultured CHN population obtained 
through artificial propagation. All cultured populations 
exhibited the less genetic diversity for gene diversity and 
nucleotide diversity compared to the wild populations. 
2.3 Genetic differentiation among populations 

Tab. 4 shows different indices of dissimilarity 
between pairs of populations, including genetic distance 
and Fy, statistics. When populations were compared 
pairwise, the highest F4 values appeared in CHN/WLJ 
pairs (F,70.2991, P«0.01). The Fx values for WLJ/WSS 
were also statistically significant (Fy=0.1414, P«0.05). 
The comparison of populations WSS and CHN, CSS 
show low F4 values, indicating that these two 
populations were not highly divergent. But interestingly, 
genetic differentiation between CHN and CSS was 
significant (Fy=0.1376, P<0.05). 
Fx values, the pairwise Nei's genetic distance value was 
low between WSS and CSS, CHN populations (0.099 


In concordance with 


and 0.298, respectively). The largest genetic distance was 
between WLJ and CHN (1.735). 


3 Discussion 


The artificial breeding program for this species 
starts for a short time and will expand accordingly. There 
is no previous genetic record of this species and no 
attempt has been made to assess the genetic status of 
both wild and cultured populations of this species. This is 
the first study to demonstrate that a mitochondrial 
marker can be used to monitor changes in the genetic 
diversity of the hard-shelled mussel (M.corcuscus) 
during domestication. 

Genetic variability is an important attribute of 
species under domestication, since those with higher 
levels of variation are most likely to present high 
additive genetic variance for production traits (Alarcon et 
al, 2004). In this study, genetic analysis of mitochondrial 
COI gene sequences revealed higher genetic variability 
for gene diversity and nucleotide diversity both in wild 
and cultured M. corcuscus than for the same sequence 
regions in other marine organism. Results from 
mitochondrial COI gene sequences in crayfish showed 
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Tab.3 mtDNA haplotypes distribution of Mytilus coruscus populations 
(GenBank accession no and number of individuals) 
Name GenBank accession No WSS CHN CSS WLJ Total 
Hapl FJ495257 1 
Hap2 FJ495258 6 7 4 1 8 
Hap3 FJ495259 2 1 1 4 
Hap4 FJ495260 1 2 1 4 
Hap5 FJ495261 1 
Hap6 FJ495262 1 
Hap7 FJ495263 1 
Hap8 FJ495264 2 5 4 1 
Hap9 FJ495265 1 
Hap10 FJ495266 1 1 
Hapll FJ495267 1 
Hap12 FJ495268 1 
Hap 13 FJ495269 1 1 
Hap14 FJ495270 1 1 
Hap 15 FJ495271 1 2 1 4 
Hap 16 FJ495272 1 
Hap17 FJ495273 1 
Hap18 FJ495274 1 
Hap19 FJ495275 1 
Hap20 FJ495276 1 
Hap21 FJ495277 1 
Hap22 FJ495278 1 
Hap23 FJ495279 1 
Hap24 FJ495280 1 1 
Hap25 FJ495281 1 1 
Hap26 FJ495282 1 
Hap27 FJ495283 1 
Hap28 FJ495284 1 
Hap29 FJ495285 1 
Hap30 FJ495286 1 





Tab.4 Pairwise estimates of Nei's genetic distance (above diagonal) and Fs 
(below diagonal) between Mytilus coruscus populations 





WSS CHN WLJ CSS 
WSS 0.09955 0.755" 0.2988 
CHN 0.0230^5 1.735" 0.954" 
WLI 0.1414" 0.2991” -0.05055 
CSS 0.04665 0.1376" -0.00875 


* P«0.05;" P«0.01; NS, not significant. 


nucleotide diversity ranging from 0.0196 to 0.4396 
(Trontelj et al, 2005). In the mud crab Scylla serrata, the 
COI sequence divergence ranges from 0.17 % to 0.46% 
and gene diversity ranges from 0.37 to 0.85 (Fratini & 
Vannini, 2002). Relatively high levels of DNA diversity 
characterized the population of the bivalve mollusc 
Congeria kusceri (haplotype diversity=0.66 in the COI 
gene) (Stepien et al, 2001). However, the mtDNA 
nucleotide and gene diversity were reduced in the 
cultured population. A slight decrease of genetic 
variability in cultured populations had been observed in 
many fish and mussel species (Lundrigan et al, 2005; 
Pampoulie et al, 2006; Wang, 2007; Kong & Li, 2007; 
Shu et al, 2008). The reduced genetic variability we 


observed in the farmed strains is probably due to a low 
number of successful breeders during the foundation 
period, similar to a recent bottleneck in terms of impact 
on genetic variability (Allendorf, 1986). The greatest 
reductions in number of haplotypes, nucleotide diversity 
in cultured CHN 
population in this study are probably caused by the use 


and gene nucleotide observed 


of small numbers of brood-stock collected from wild 
The in the cultured CSS 
population are probably contributed to collection of 


populations. reductions 
natural spat from small numbers of parents and samples 
examined. A more precise assessment of the genetic 
variability in cultured populations can be made, were the 
magnitude of the genetic variation of wild populations 
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made available. This is almost a requirement when the 
number of generations of the cultured populations is 
small and there has not been enough time for the 
variation to be reduced to a detectable level. 

The comparison of the Ft values (Tab. 4) between 
CSS, CHN and WSS populations showed that the values 
were not significantly different. This indicates that no 
significant genetic differentiation was detected between 
CHN, CSS and WSS populations. Similar findings were 
reported for fish (Yang et al, 2008), marine mussel (Jiang 
et al, 2007; He et al, 2008). The slow genetic 
differentiation might be due to persistent gene flow 
caused by cultured juvenile practice in open sea, close to 
the wild population. But significant genetic 
differentiation observed between WLJ and WSS, CHN 
populations, was probably associated with geographic 
factors. This regional variation could be an important 
source of diversity for both genotypic and phenotypic 
traits to be selected in accordance with aquaculture goals 
(Lundrigan et al, 2005). 
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上 海 海洋 大 学 水 产 动物 种 质 资源 与 创新 团队 简介 


水 产 动 物种 质 资源 与 创新 团队 是 上 海 海洋 大 学 水 产 养殖 国家 重点 学 科 的 骨干 研究 力量 ， 也 是 省 部 共 建 
水 产 种 质 资 源 发 据 与 利用 教育 部 重点 实验 室 的 基本 研究 队伍 。 


上 海 海洋 大 学 是 国内 最 早 开展 水 产 动物 种 质 资 源 与 创新 研究 的 单位 ， 世 界 水 产 养殖 终身 成 就 奖 获 得 者 
一 一 李 思 发 教授 80 年 代 初 回国 后 创建 了 这 个 研究 领域 ， 至 今 已 有 近 30 年 的 历史 ， 在 国内 外 产生 了 广泛 的 
学 术 影 响 。 李 思 发 教授 领导 的 课题 组 历时 20 年 经 过 6 代 系 统 选 育 , 获得 了 生长 速度 比 原 种 快 30%、 体 形 好 、 
遗传 性 状 稳 定 的 团 头 芒 新 品种 “浦江 一 号 ”， 它 是 世界 上 第 一 个 经 选 育 的 草食 性 鱼 类 新 品种 ， 也 是 上 海 市 
育成 的 第 一 个 水 产 动 物 良种 ， 该 项 目 获 2002 年 度 上 海 市 科技 进步 一 等 奖 、2004 年 度 国 家 级 科技 进步 二 等 
奖 。 他 领导 的 课题 组 还 于 2006 年 初 选 育 出 新 吉 富 罗 非 鱼 新 品种 ， 以 该 新 品种 为 核心 的 “从 吉 富 到 新 吉 一 
一 尼 罗 罗 非 鱼 种 质 创新 与 应 用 ”项 目 荣获 2007 年 度 上 海 市 科技 进步 一 等 奖 。 

该 团队 年 轻 学 术 带 头 人 李 家 乐 教授 是 李 思 发 教授 的 学 生 ， 他 领导 的 课题 组 从 1998 年 开始 ， 对 我 国 具 
有 悠久 养殖 历史 、 产 量 占 世 界 产 量 95% 以 上 的 淡水 珍珠 蚌 ， 开 展 了 系统 的 种 质 评价 与 第 选 ， 获 得 了 三 角 帆 
蚌 性 阳 湖 新 品系 ， 并 率先 将 三 角 帆 蚌 和 从 日 本 引进 的 池 蝶 是 进行 杂交 ， 获 得 了 新 品种 “康乐 蚌 ”。 “淡水 
珍珠 蚌 新 品种 选 育 和 养殖 关键 技术 ”项 目 获 2008 年 度 上 海 市 科技 进步 一 等 奖 。 

目前 ， 该 团队 拥有 教师 16 名 ， 其 中 教授 6 人 ， 副 教授 5 人 ， 是 一 支 以 中 青年 和 具有 海外 留学 经 历 的 
教师 为 主 的 创新 型 团队 。 该 团队 拥有 研究 生 75 名 ， 其 中 博士 研究 生 12 人 。 

最 近 五 年 来 , 该 团队 获得 3000 多 万 科研 经 费 ， 其 中 重要 的 科研 项 目 有 : “973” 前 期 研究 专项 、“863”、 
科技 支撑 、 农 业 产业 体系 、 国 家 自然 科学 基金 重点 等 。 获 国家 科技 进步 二 等 奖 1 项 ， 上 海 市 科技 进步 一 等 
X 3 项 ， 省 部 级 二 等 2 项 、 三 等 奖 10 余 项 ; 在 公开 的 学 术 刊 物 上 共 发 表 论文 253 篇 ， 其 中 SCL. EHR 
26 篇 ; 主编 或 参 编 专 著 5 部 、 教 材 5 部 ， 其 他 编著 6 部 ; 申请 专利 12 项， 其 中 获得 3 项 。 
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